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, We study the (local) propagation of plane waves in a relativistic, non-dissipative, two-fluid system, 

' allowing for a relative velocity in the "background" configuration. The main aim is to analyze 

^ ' relativistic two-stream instability. This instability requires a relative flow — either across an 

interface or when two or more fluids interpenetrate — and can be triggered, for example, when one- 
^N) ■ dimensional plane-waves appear to be left-moving with respect to one fluid, but right-moving with 

' respect to another. The dispersion relation of the two-fluid system is studied for different two-fluid 

equations of state: (i) the "free" (where there is no direct coupling between the fluid densities), (ii) 
coupled, and (iii) entrained (where the fluid momenta are linear combinations of the velocities) cases 
are considered in a frame-independent fashion [eg. no restriction to the rest-frame of either fluid). 
As a by-product of our analysis we determine the necessary conditions for a two-fluid system to be 
^H^. causal and absolutely stable and establish a new constraint on the entrainment. 
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I. INTRODUCTION 



Newtonian physics is replete with examples of multi-fluid systems, such as diffusion, ion flow, superfluid Helium, and 
plasma discharge from the Sun. In fact, large characteristic scattering times between different components is more 
■ the norm than the exception. This leads to physical situations where the various components can move independently 
• ■ of each other, be it across an interface or through interpenetration. In this context, even heat conduction in systems 
: where all the matter flows together is a two-fluid problem, i.e. there is a heat flux in addition to the matter flux. 
Qv^ ' Perhaps not as widely appreciated is that the relativistic regime has its own set of multi-fluid scenarios: neutrino 
, streaming during supernovae, superfluid neutrons and superconducting protons in neutron stars, and heat flow in a 
1!^ ' cosmological setting, to name but a few. 

A key issue is that relativistic fluids must be causal, meaning that sound speeds, say, must be less than that of light. 
For fluids, there are two entry points for causality: the microscopic where particle-particle interactions are tracked 
5— i ■ and the macroscopic where fluid elements (large enough to contain many particles, but small enough to be point-like 
with respect to the total system) are monitored. Presumably, a fully relativistic treatment at the microscopic level 
would lead to a set of fluid coefficients (describing the equation of state, dissipation, etc.) that would already behave 
appropriately at the macroscopic level. However, there is a practical problem: Equation of state determinations 
are notoriously difficult. This makes a general analysis of relativistic fluid dynamics prohibitive, if not impossible. 
Fortunately, one can make progress by imposing causality "from above" and absolute stability {i.e. real sound speeds) 
"from below" to constrain the fluid coefficients. 

In this paper, we will do this by analyzing the local propagation of plane waves on a given (arbitrary) background 
spacetime. Compared to the standard single-fluid analysis, we have more fluid degrees of freedom and need to 
allow for relative flows between the various fluids. This is an essential requirement for two-stream instability. Such 
instabilities arc known to exist for a variety of configurations. For shearing motion at an interface, it is an example of 
Kclvin-Helmholtz instability. However, as far as we are aware, generic two-stream instability has not been discussed 
previously in relativity. 

The two-stream instability has been well-documented for plasmas (where it is known as the "Farley-Buneman" 
instability [l|, see 0] for a text-book discussion). It has also been suggested as the mechanism behind star 
formation when two galaxies (whose angular velocities are more or less anti-aligned) merge In the general 

relativistic context, Chandrasekhar, Friedman, and Schutz (CFS) 0, Q have demonstrated that oscillation modes in 
rotating, perfect fluids can become two-stream unstable due to the emission of gravitational radiation. Here, the two 
"fluids" are the rotating mass, and the asymptotically flat spacetime in which the fluid is embedded. Most recently, a 
two-stream instability for supcrfluids has been proposed, with a natural extension to a mixture of superfluid neutrons 
and superconducting protons in neutron star cores 0, Q . Very recent results suggest that this instability may act as 
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a trigger mechanism for the enigmatic spin gUtches that have been observed in a number of radio pulsars [9(] (see also 
Q for the first suggestion of a link between glitches and two-stream instability). 

In what follows we will not restrict the discussion to any specific physical system. Consequently, the analysis 
will be somewhat abstract. This strategy can work because the two essential requirements for triggering two- 
stream instability is a relative flow between two fluids and a generic interaction between them. This freedom to 
remain abstract illustrates the general robustness of the instability and its presence in a diverse collection of systems. 
Essentially, if the relative velocity is large enough that a wave moves in one direction with respect to the rest-frame 
of one fluid, yet the opposite direction with respect to the other fluid's rest-frame, then the energy of the wave will 
be "negative" in one of the rest-frames and therefore unbounded from below. Our main aim is to show that a causal 
and absolutely stable system of two relativistic fluids can undergo two-stream instability for a range of relative speeds 
and equation of state parameter values. 

The presentation of the results is organized as follows: Section |TT] recalls the multi-fluid formalism, and sets the 
stage for a plane-wave analysis of the system. Section IIIII considers sound waves for a single fluid. The results 
are not new, but they help establish basic techniques that carry over to the more complicated two-fluid calculations 
discussed in Section llVl The following sub-sections consider three variations on the two-fluid equation of state: (i) the 
"free" (where there is no direct coupling between the fluid densities), (ii) coupled, and (iii) entrained cases. Finally, 
m Section |Vl we make our concluding remarks. Spacetimc indices are denoted by the first letters of the roman 
alphabet (a, 6, c), constituent indices by the last (x, y, z), and we adopt "MTW" (Misner, Thornc, and Wheeler [To| ) 
conventions for the metric signature. 



II. THE MULTI-FLUID FORMALISM 

We will use the approach to multi-fluid systems that was developed originally by Carter [ll| (see [l^l for a recent 
review). In this description, the main variables are the various fluxes (for particles and/or entropy), to be denoted 
and the equations of motion follow from a suitably defined "master" function (i.e. Lagrangian or equation of state) 
A. In the single fluid case, —A is equal to the rest frame energy density p . We have here introduced the convention 
of attaching a constituent index x to each variable. This index is redundant for a single fluid, but necessary when 
there are multiple fluids. The master function varies only with the fluxes. If the fluids are locally isotropic [i.e. no 
preferred direction), as they should be in the absence of anything else (such as an elastic solid), it is clear that A must 
be a function of only the various scalars that can be formed from inner products of the fluxes. 



We will focus on the case of two fluids (see [l2| for a complete description), even though most of the equations in 
the general discussion will carry enough constituent indices (x, y, etc.) to be valid for any number of fluids. In the 
case of two components, the master function depends on two distinct particle fluxes n'^ and and has the functional 
dependence [1^ 

A = A{nl,nl,nly) , (1) 

where n^y = —gabn'^n^- Note that constituent indices are not summed over when repeated. As a matter of 
convenience repeated indices are written only once; that is, we write (which is the squared particle number density 
of the x"^-fluid) for n^^ and so on. 

The equations of motion become most transparent when expressed in terms of the momentum which is canonically 
conjugate to n^: 



where 



l,l = B-nl+A-^nl, (2) 

^ -^-^^--T^i- (3) 

Note that we have simplified the notation by not indicating explicitly which variables are fixed when partial derivatives 
are taken. (The functional dependence of the master function is clear from ([1]).) From ^ we see that the momentum 
/i^ is not simply proportional to its canonical conjugate n°, but is rather a linear combination of all the fluxes. This 
is a result of the so-called entrainment effect (see for an example in superfluid Helium mixtures, [ij, [l^, [l^ for 
relativistic, nuclear matter, or p^j for a treatment of entropy /matter entrainment and its importance for heat flow). 
It is convenient at this point to introduce a shorthand notation for derivatives of these coefficients; namely, 

^ ^ (2n n ^ 
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SI, 



•^,xy 



dnly 



For the same reason we define the "speed-of-sound" of the x"^-fluid as 

2 d log 



9 log TT-x 



Finally, we introduce the "perp" operator 



\ xb r & I X 6 
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(4) 



(5) 



(6) 



which can be used to construct, say, vectors that are orthogonal to 

We have chosen the fluxes n° as our primary fields. However, there is no reason why the momenta could not 
be similarly adopted [2^ . This implies that the mapping from one set of fields to the other must have an inverse. 
That is, we see from ([2]) that 



(7) 
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and thus B^^B^ 
sound modes. 



(A^^) 7^ 0. This will be a useful constraint later, when we discuss the impact of entrainment on 



For the purely variational case (i.e. no dissipation, imposed constraints, etc.), the individual constituents are 
|12l . [l8| , so that we have for each component 



conserved 



The remaining equations of motion take the form 



where the vorticity tensor oj^^ is given by 



Va< = 



2V 



afJ-b] 



(8) 
(9) 
(10) 



As discussed, eg. , in [12l |. in the single-fluid case these equations contain the same information as the standard set 
obtained from the vanishing of the covariant divergence of the stress-energy-momentum tensor. Eq. © illustrates 
the geometrical significance of the Euler equation as an integrability condition on the vorticity; i.e. that the particle 
fiux nowhere pierces the surfaces defined by the two-form w^f,. 

Below we will be analyzing plane- wave propagation on backgrounds such that = 0, the various background 
quantities are taken to be constant, and there is a relative flow between the fluids. This implies a linearization of the 
equations of motion; i.e. 



Va6K = 



?i =0 



(11) 



Because there are several fluids, the variation is signiflcantly more complicated than, say, in the case of the 
standard perfect fluid. In addition to individual bulk effects for each fluid, there can also be cross-constituent effects 
due to coupling between the fluids. We also have to consider the entrainment. 
Following |12l ]. we can isolate the various effects in the variation and write 



s^l: = (b:, + a:,) sni + {x:i + az) sr 



where the bulk effects arc captured by 



the cross-constituent coupling through 



Bib 



B'-il-lb 



Cx'"a'"b) 



(12) 
(13) 



-CccVB^uW 



(14) 
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and the entrainmcnt via the terms 




In these expressions the flux n'^ has been decomposed as = n^u'^, where u° is the unit (u^u'^ = four- velocity 
of the X- fluid elements. 



III. SINGLE-FLUID SOUND WAVES 



To set the stage for the general analysis it is useful to first consider the nature of sound waves in a single, perfect 
fluid. To do this, we perform a local analysis of linear perturbations of the fluid (keeping the metric fixed) on a 
generic background. In particular, plane-wave propagation corresponds to the Ansatz 

Sn^ = exp(ifcf,x'') , (17) 

where the amplitude A° and wave four- vector have vanishing covariant derivatives. Recall that we assume all 
unperturbed quantities to be similarly constant. In particular, the background vorticity simply vanishes. From p2|) 
above we have 

Sf^l = B^M ■ (18) 

Note that we have continued to use a constituent index, even though we are dealing with a single fluid. This allows 
for some economy of presentation since many of the formulas will apply later, except that the index x will then range 
over two fluids. Of course, waves in a system are such that the constant wave vector ka is the same for all the fluids. 
Hence, it does not carry a constituent index. 

For convenience we will work in the material frame associated with the fluid. This means that ka and A'^ will each 
be written as two pieces by utilizing the "perp" operator introduced in ([6]): 

< + k:) , (19) 
where cr^ and the wave vector fc^ (with magnitude fcx) are 

k-x_(^x = ^kaU^ , k^ =-Lxa = ^x^a ■ (^0) 

Similarly, we can decompose the wave amplitude as 

Al=A^^ul + Al^, (21) 

where 

Al = -ulAl , Al^=l.l,Ai. (22) 

Note that CTx and are measured by an observer moving with the fluid. It will be obvious from the dispersion relation 
constructed below that Cx measures the phase velocity of the waves as seen in the fluid frame. Furthermore, it is 
easy to show that evaluating Cx in a frame moving relative to the fluid leads to the standard Lorentz transformation 
of velocities. 

With these preliminaries, the perturbation equations (jlip reduce to 

= kaAl , (23) 
= Kk^aBl^^Al . (24) 

The first of these relations shows that the waves are transverse in the spacetime sense. The dispersion relation can 
be easily obtained by contracting the second equation with fc''. (In the more complicated two-fiuid analysis below, we 
will in general have to consider the vanishing of a 4 x 4 determinant.) Assuming that n^A'^ ^ 0, and after several 
steps of algebra, the dispersion relation reduces to 



(25) 
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In our homogeneous plane- wave setting it is clear that the group and phase velocities coincide so that we can introduce 
the speed of sound in the standard way as = a^. 

To see that this is equivalent to the usual single-fluid result (as in, say, [3), it will suffice to introduce the pressure 
and recall that p = —A in the single fluid case. From the definition above for B^, cf. we see that 

dp = B^'n^dn^ . (26) 

Moreover, the pressure p is defined by the standard thermodynamic relation, such that 

p = -p + n^^ = -p + B^'nl . (27) 



This implies 



d log rix 
and we have 

dp d log B^ 



'dnx 



dP=(l + l!^MP, (28) 



4 , (29) 



dp d log rix 



where dp/dp is the usual form of the sound speed squared. 

In order to pave the way for the more complicated multi-fluid case to be discussed below it is useful to examine the 
properties of the various vectors we have introduced. Starting with the wave vector fc" we see that 

kak'^^kKl-cl) . (30) 

Thus, for causal wave propagation (c^ < 1), fc" is spacelike. For the wave amplitude we find that, when the force 
equation is evaluated in terms of the solution to the dispersion relation, 

A»^=a.Aj^fc^ (31) 

The waves are therefore longitudinal in the normal, three-dimensional sense. On the other hand, the transverse 
nature (j23p of the waves in the four-dimcnsional sense implies that 



Al = Al^ (1 - c-') (32) 

so that A^ is timelike (and we can choose it to be future pointing) for causal waves. Note that since the flux is 
also timelike this implies that n^A^ < and thus the degenerate case of ([25]) is ruled out by causality. 

Before we conclude this section, it is useful to consider what constraints (O imposes on the equation of state. First 
of all, the causality requirement leads immediately to 

t^ + l<l. (33) 

dlog Tlx 

In addition, we must have in order to avoid absolute instabilities (complex wave speeds). This implies 

d log B^ 



9 log Tlx 

Combining the two results we see that we must have 



-fl > . (34) 



-1<§;^^<0. (35) 
d log rix 



In the next section we will extend this type of analysis to the two-fluid model. 
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IV. SOUND WAVES FOR THE GENERAL TWO-FLUID SYSTEM 



We want to work out the dispersion relation for wave propagation in a two-fluid system, using (|17[) as the starting 
point. An important addition to the problem of plane waves is a new "parameter", the relative flow between the 
two fluids. We will represent this flow by the relative velocity of the y'^-fluid with respect to the frame of the 



^th 



fluid: 



where v-^y represents the magnitude of the relative flow and 



(36) 



7xy 



c y 

7yx — ^'^viU^ 



(37) 



This leads to 



h = 7xy [U 



(38) 



The mode speed CTx and wave (three-) vector can be deflned as before. The insertion of 7xy into (jSG]) makes the 
x*''-fluid's proper time the standard for setting velocities. Since the dispersion relation below is a scalar equation, 



we will have in several places the inner product v^yk^ (where v^y 
the angle O^y between the two vectors: 



v^y/vxy)- It is useful to write this in terms of 



''xy'-a 



cos( 



xy 



(39) 



An important subtlety must be recognized, however: The three quantities CTx, k^, and v^^ are what would be 
measured by an observer flowing with the x"^-fluid. We could have equally as well chosen the material frame 
attached to the other fluid (or some other observer). As one might expect, there are well-defined transformations 
between the two descriptions (which will be needed later). The relative flow of the x"^-fluid with respect to the 



y"^-fluid frame is related to via 



V. 



yx — 7xy (''^xy^x + ''-'xy) 



where we have used v- 



yx 



^xy 



It is also useful to note that 

— ^V^y (Vxy 

«y = -«xy'«x 



/xy "yx; 
■7"^W° ) 



Because ka is not attached to cither fluid frame, we must have 



ka = ky (ayul + klj = fcx (cr^ul + k^^ 



By contracting each four- velocity in (|4ip with the wave- vector ka, we obtain a matrix equation for [k^ ky]"^ 



VxyCTx - COS^xy 


-7x/cos6'yx 




kx 




■ 0" 


-7x/cos6ixy 


VxyCTy — COS Syx 




ky 








(40) 

(41) 

(42) 
namely, 
(43) 



Obviously the determinant of the 2x2 matrix must vanish. This leads to 



cos 6'v 



Wxy cos 6'xy 



cos ( 



xy 



(44) 



It is not difficult to show that if cr^ < 1 then ay < 1. This is natural given that causality is a frame-independent 
requirement. 

Meanwhile, the equation of flux conservation is the same as ([23|) (except x ranges over two values). The conservation 
of vorticity equations become 



^ab^y 



7>-xy .f, 
^ab^y ' 

j.-yx /if- 



(45) 
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where the "dispersion" tensors are 



Kl = <(fc[,A'^f, + fc[e-4^fJ . (46) 



Note that ivT^j, and arc obtained via the interchange of x y in equation (|46 
In order to solve (|45p. we obviously need the four inverses 



K-^K^,^5\ , K-^K^^ = S\, (47) 
to exist — i.e. the determinants of K^i^ and A'^^ do not vanish — so that we can write 

= [k^'KI^ - K^^K,) Ai ^ MabAi . (48) 
The only way to get a non-trivial solution is to have a ka such that 

e-^---^''^e''''^-''''''Ma,b,Ma,b,Ma,b,Ma,b, = . (49) 

Written out in full ((49|) is a quite busy expression. This should come as no surprise since the two-fluid problem is 
significantly more complicated than a single fluid having bulk contributions coming from S^^, and the cross-coupling 
from X^^ (in the case of two, co- moving constituents [l^l)- For two-fluid systems, the two constituents move 
independently and there are the additional contributions A^f^ and A^ coming from entrainment. 

In order to simplify the problem, it is convenient to isolate further the different contributions that appear in the 
dispersion matrices. The bulk contribution in K^f^ can be reduced to 

Kb = Kk[cBa]b 

= -l-B^n^{k^a^±:,+clk^,ul) , (50) 



while its entrainment piece becomes 



n. 

Meanwhile, the cross-coupling in K^^ is simply given by 

^ab — ^^^lc^a]b 
1 



<hc^l]b 

"^AJ^ {k^a^v^y - kl) {ut + <")} . (51) 



= --C,,7xyn.V6^^'^« + <^) , (52) 



but the entrainment has significantly more presence: 



A y + ^ ^.xy + 7xy ( 7xy + •^.xy 

y 



+7xyfcxax (7xy^ef,y + Aj)j ^^K^j . (53) 

A multi-fluid system must have non-zero bulk properties (unless a fluid vanishes completely). The other terms 
can be absent, depending on the equation of state. In what follows we will systematically increase the complexity by 
considering in turn the different components of two-ffuid physics. 
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A. Dispersion Relation: Free Case 

Let us first consider the case of two completely uncoupled fluids. Then we have only non-zero. In lieu of (|49p . 

it is easier to get the dispersion relation by contracting the free indices in (I45p with ka- This results in the simple 
2x2 matrix problem 



{al ^4) 








■ 0" 


By {al - 4) _ 












(54) 



and the dispersion relation is simply 

(a^ - 4) i'^'y - 4) = . (55) 

By construction cr^ is the squared phase (three-) velocity as measured in the x'^-fluid frame (similarly for <Ty). 

The outcome of this analysis is that the dispersion relation (|55p allows four non-trivial solutions consisting of the 
roots of 

al = 4 (56) 

for any x. These roots correspond to ±Cx evaluated in the x*'^-fluid frame and are just Lorentz transformed if 
evaluated in another frame. Thus, as expected in the case of zero coupling between the fluids, the quantity Cx can 
be interpreted as the sound velocity of the x^'^-fluid as measured in its own (background) rest frame. This is the 
obvious generalisation of the single fluid result. It follows that if Cx is subluminal in its own rest-frame it is so in all 
other frames as well. Also, since absolute instability should be evaluated at zero relative velocity it is clear that the 
constraint < Cx < 1 remains as a condition for the master function. 

The main conclusion from this discussion is that the constraints on the equation of state are easily generalized 
to the uncoupled two-fluid problem. We must thus require that the equation of state satisfy ([35]) for both fluids 
in order for the system to be absolutely stable and give rise to causal wave propagation. There are no dynamical 
instabilities present in this case. In what follows, we caution that the Cx can be understood as "sound" speeds only in 
this completely free case. When fluid couplings are operative, the phase velocities will no longer simply equal these 
free sound speeds. But in order to make progress, we will continue to impose psp throughout. 



B. Dispersion Relation: Cross-constituent Coupling 

We now allow for to be non-zero. From ([5^ . we see that this will introduce Ccc in addition to B^ and 4 
Using the same contractions with ka as in the free case, we again get a 2 x 2 matrix problem; i. e. 



B''{<jI-4) -VB^ByCcc 








' 0' 


-VB-Byc,, By (a2 _ 4) \ 












(57) 



The dispersion relation is now 

{al - 4) {al - 4) = Cl . (58) 

In order for to be less than zero, the equation of state would have to allow either of the B^ to be negative. For 
ordinary matter, or entropy, this is not generally the case. Hence, we do not consider this possibility here. We also 
need to point out that (|33]) must be used to get a dispersion relation solely in terms of CTx- Clearly, the cross-coupled 
case is more complicated than the free problem. However, it is also much more interesting and relevant. As we will 
soon see, the richer phenomenology allows for two-stream instability. 
Addressing first the question of absolute stability we set v^y = to find 




In order to avoid complex a'^, the discriminant of (j59p must be positive, which is evident for > 0. Recall that 
absolute stability means CTx > 0. Clearly, it is sufficient to require that 



{4 + 4)">i<^l-4)' + ^i 



(60) 
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The second term of ([59)1 will always be less than the first if 



In other words, since we expect to have B^B"^ > 0, absolute stability constrains the equation of state to satisfy 



dlogB"" dlogB^ 
dloguy 9 log Tlx 



< 1 



aiogs^ 

d log Tlx 



1 



d\ogB^ 

d log Uy 



(61) 



(62) 



The causality constraint requires cr^ < 1- For C^^ > and satisfying (|6T|) . we need only make the "- 
([55)) causal, since the "— " solution is always smaller. We find that causality is ensured if 

Cl < (1 - ci) (1 - 4) , 

which, in terms of the equation of state, translates into 

dlogB"" dlogBy 



solution in 



(63) 



dlogB"" dlogBy 
91og% d log Tlx ~ 9 log fix 91ogriy 



< 



(64) 



Note that both and restrict C^^ from above. In general, we can show that if + Cy < 1 then any absolutely 
stable equation of state is also causal. Conversely, if + > 1, any causal equation of state is absolutely stable. 

Given an equation of state that is causal and absolutely stable, we can now determine if a dynamical two-stream 
instability is present by solving the dispersion relation (|58p for some relative flow (as parameterized by Wxy). Writing 
down the general solution is not difficult, but it is instructive to first focus on the slow velocity limit. Assuming that 

, the dispersion relation (|58p becomes 

(ax-«xyCOS0xy)'-d =C2^ . (65) 



and (Tx are both much smaller than unity (the speed of light) in 



Introducing new variables 



Ox 
Cv 



y = 



l^xy COS 0xy 



(66) 



we get 



-62 



(x-y) -1 



= 1 



(67) 



As one might have expected, the problem is now identical to the Newtonian plane-wave problem discussed by Ander- 
sson. Comer, and Prix Q. Hence, we can learn from their results. They demonstrate that a two-stream instability 
may operate above a critical relative flow. Their particular example corresponds to = 0.0249 and = 0.0379. 
For this case they flnd an instability in the range 0.6 < y < 1.5. This means that the system becomes unstable for 
Cxy > 0.6. This flow is clearly sub-luminal as long as Cx < 1, but one may suspect that the linear approximation 
that we have used is not very accurate. Still, this is a useful first demonstration that the two-stream instability will 
operate also in relativistic systems. 

Before we turn our attention to the full relativistic case it is useful to check if the particular example used by 
Andersson et al. Q obeys the causality and absolute stability criteria derived above. First we note that, due to the 
presence of a velocity scale given by the speed of light, in relativity we cannot completely scale out the velocities. 
Thus the relativistic analysis of stability will in general contain an extra parameter compared to the Newtonian case. 
Here we shall take that parameter to be Cy which, without loss of generality, can be taken to be larger than Cx- Using 
these parameters the absolute stability criterion ([6T|) becomes just 



a 



(68) 



which is satisfied in the model discussed above. The causality condition ([35|) enters only indirectly as we have been 
able to re-scale in terms of Cy and thereby get dimensionless variables. We conclude that the Newtonian model of 
Andersson et al. is reasonable also from this perspective as long as Cy is not very close to the speed of light. 

We now turn to the relativistic dispersion relation (|59[) . Written as an equation for it constitutes a non-trivial 
quartic. If we use the same re-scalings as in (|66|). then ([59]) becomes 



"fxy 



CyX^ 



) +cl{x~y) 



1 



1 . 



(69) 
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FIG. 1: Plots of the real (top) and imaginary (bottom) parts of the mode frequencies as functions of v^y, for Cy = 0.5, 
— 1.0, and Sxy = 0. The merger of two frequencies, and subsequent non-zero imaginary values, signal the presence of a 
two-stream instability. 



Some immediate insight is obtained by considering the ultra-relativistic Umit for the background flow, where v^y 1 
or 7~y^ — ^ 0, and the hmit where the wave vector becomes perpendicular to the background flow, i.e. ^xy ^ 7r/2. In 
both hmits the two-stream instability ceases to operate. 
In the ultra-relativisitic limit, the wave-speed tends to 

[ cos6'xy (double root) 

If the propagation is to remain causal, we must have C^^ — s- as Cy — > 1. Also there is no two-stream instability 
since C^^ > 0. This might seem surprising, since a two-stream instability requires a "window" of background flows 
for modes to appear, say, left-moving in one frame but right-moving in the other. But as — > 1 the relative flow 
is at its maximum, and yet the instability window is closed. In fact, this behaviour was seen already by Andersson 
et al. [1] in the Newtonian regime. From we also see that y — *■ as 6'xy — > 7r/2. This turns into a quadratic 
for , and one finds that the discriminant is positive for the range of values for C^^ that yield absolute stability and 
causality. Obviously, y is the effective "window" of the background flow and it is completely closed for O^y = 7r/2. 

Although the general solution to is readily availiable, it is quite complicated and offers very little additional 
insight. Instead of writing it down we will tackle the problem numerically. The parameter values are restricted 
to those that maintain absolute stability and causality. Figure [1] provides plots of the real and imaginary parts for 
the four solutions to ((55]) in the aligned case. The solutions for CTx are taken to be functions of the relative flow 
parameter y and the coupling C^^, with Cy = 0.5, = 1, and ^xy = 0. Non-zero values for ImcTx in the figures 
indicate the presence of an unstable mode. The appearance of unstable modes is reflected in the real parts wherever 
two frequencies merge. This behaviour is typical for this kind of dynamical instability. The results for misaligned 
flows, with 6'xy 7^ are very similar to those shown in Figure [T] As 6'xy increases the [y, C^^ region of instability 
moves towards higher relative velocities, eventually leading to regions that stretch essentially all the way to v^y = 1. 
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C. Dispersion Relation: Aligned or Anti-aligned Flows 



Now that we have estabhshed the presence of the two-stream instabiUty for arbitrary background flows, we will 
consider the more restricted case of aligned or anti-aligned background flows. This will simplify the dispersion relation 
so that analytical insight can be more easily acquired. It also reduces the parameter space, thus allowing a more 
focused numerical analysis. 

By aligned or anti-aligned flow we mean that the wave propagation is aligned (or anti-aligned) with the relative 
velocity of the two fluids; specifically, 0xy = (0, tt) so that 

K^e^yVa^ , Cxy = COS^xy = ±1 . (71) 

This leads naturally to the statement that the wave vector is a linear combination of the background flows: 

k 

ka = — [(WxyO-x " Cxy) < + ^xyTxy ^^a] 
"Wxy 

k 

= — [(WxyO-y - eyx) < + eyx7x/Ua] ■ (72) 
"xy 

Equating coefficients in ([TS]) leads to 

ka = {kyEyxU^ -\- fcxExy^ia) ■ (73) 

7xyWxy 

Because of ([75)) . our original four-dimcnsional linear algebra problem for {u^A'^, it^Ay , u-^A^, u^Ay)"^ has been reduced 
to a two-dimensional one for {u^A'^,u^Ay)'^. 



D. Dispersion Relation: Role of Entrainment 



Up to this point we have introduced three equation of state parameters (cx, Cy, and Ccc) that are obtained as second 
derivatives of the master function. When entrainment is included in the model we see from Q that we need two 
additional variables to describe the general case. Given that we have established the two-stream instability for general 
cross-constituent coupling and arbitrary background flow, the main reason for discussing the role of the entrainment 
is to highlight the basic feature that the instability can be triggered by a variety of interactions. We will simplify 
the entrainment case by assuming that the relative velocity v!^y is much smaller than the speed of light and that the 
flows are aligned (in the sense of the previous section). This docs not mean, however, that the individual flows u° 
have to be similarly restricted. Neither do the sound and wave speeds Cx and ctx have to be small. 

If we keep the relative flow to 0{v^ ), the master function can be approximated as pol. [2l| 



A = Ao(7ix>7iy) + Xi{nl,ny) (^nly - 



lly 



(74) 



which immediately implies 



B"" = -2 



9Ao ^ dXi f 



6^5 



9 ^ 



(75) 



and 



/ixy 
•^,xy 



We shall make one further simplifying approximation, which is to take Ai to be constant so that 







(76) 
(77) 



This leaves us with the single entrainment parameter Ai. 
With these approximations wc find 



eyxfcy 



-CccVB^-A^y'i^{al 

2\ 



1) 





■ <A^. ' 




' 0' 




_ <A- 








(78) 
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From dSOl) and (HI) wc have 



l-ai 



(79) 



and thus the dispersion relation becomes 



(1 - al) (1 - a2) 



(80) 



where 



(81) 



and because the inverse of ^ must exist, \Cent\ ^ 1- Note that, since <t^ ^ 1 as 1, Cent does not affect the 

ultra-relativistic Umit. Hence, equation ([70)1 [with cos(6'xy) = ±1] is stiU vaUd in this case. 

We stress that, unUke the simpler cases, (|80|) does not hold for arbitrary propagation direction with respect to the 
relative velocity. This is important qualitatively and quantitatively, since it mirrors the fact that entrainment enters 
the master function in a fundamentally different way: At first-order in the relative velocity squared. In the dispersion 
relation, however, entrainment contributes even in the limit of zero relative velocity, because there are still two sets of 
interacting sound waves. In fact, we will now follow the earlier analysis of causality and absolute stability by taking 
this limit. 



From (l44l) we see cr^ 



CTy so that 



[Ccc+Cent 







(82) 



It is particularly instructive to consider the entrainment alone, i.e. set Ccc = 0. The corresponding dispersion relation 
is a quadratic in a^, and has the solutions 



+ 

'--X ^ 



+ 



4C2 



1 



1 



1/2 



2(1-CLJ 



(83) 



For Cxy < 1, the discriminant is obviously positive and hence the cr^ are real. In order to analyze absolute stability 
and causality we need to consider the ranges < Ce„t < 1 and 1 < Cg„j separately. We will look at Cg„( > 1 first. 
The first step is to re-write ([55]) so that the denominator is positive: 



2 

ent 



l-ci 



1/2 



2 (Cent 



1) 



(84) 



Since we are imposing c^y < 1, and Ce„t > Ij the terms outside the square root in the numerator are positive. 

" . - - . - - As for the "-" 



Hence, Cl^^ > 1 



Therefore, the solution is absolutely stable. But, we can also show that it cannot be causal 
solution, we can easily show that it is absolutely stable only if C^^j < 1, which cannot be satisfied, 
does not lead to both absolute stability and causality, and is therefore ruled out. 

The range < Ce„t < 1 is a different story, because the terms outside the radical in the numerator of ([SS]) are not 
of any definite sign. This affects the absolute stability analysis more than the determination of causality. In fact, 
the causality analysis is sufficiently straightforward that we will simply state that this requirement is satisfied for this 
range of Cg„(. In order to assess absolute stability, it is useful to introduce 



1/2 



(85) 



This allows the numerator of (j83[) to be rewritten in such a way that the absolute stability condition becomes 



- [±2 (1 - cl) (1 - 4)] + [d (1 - cl) +4(1- d)] [2 - (4 + 4)] < , (86) 



where the "±" corresponds to that of 



"y I • y 
The final step is to factorize 



{-±[2-(4 + 4)]}{-t[4(i 



4(1 



and thereby obtain 

ci)]}<o, 



(87) 



13 




where if the "+" is taken from the first factor then the "— " must be taken in the second, and vice versa. In either 
case, the factor that has the "+" is positive definite, and so the other factor must be less than zero. When the first 
' the incquaUty leads to 



factor takes the 
< rV^ 

^ent — ^yi^y 

To summarize 
this range unphysical. Meanwhile, for < Cg„j 
earlier definitions, this translates into 



2 

ent 



< 1. The other choice leads to the more restrictive condition of 



we have shown that when Cg„( 



> 1 

< 



there is cither no absolute stability or causality, which makes 
the system is causal and absolutely stable. In terms of our 



d log Tlx 



d log By 

d log Uy 



as a constraint on the master function, when the relative speed between the two fluids is sufficiently small. 

Finally, we turn to a numerical/graphical analysis for exposing the two-stream instability due to entrainment 
coupling. As in the cross-constituent coupling case, we use the re-scalings of (|66p . except that Cg„j replaces C^^ in 
a^. Equation (|80|) becomes 



1 



(Cy^' - 1) 



Cy {x - y) 



1 



(89) 



The parameter values arc restricted to those that maintain absolute stability and causality. Fig. [2] plots the real and 
imaginary parts of the four solutions to (|89|1 . As before, the solutions for ctx arc taken to be functions of the relative 
flow parameter y and the coupling Cg„j, with Cy = 0.5 and 6^ = 1.0. Fig. [2] is not so dissimilar from what wc find for 
cross-constituent coupling, thus highlighting that the instability is not sensitive to the type of coupling between the 
two fluids. 
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V. CLOSING REMARKS 

There are several examples of relativistic systems that require multi-fluid dynamics for qualitative understanding 
and quantitative accuracy. In such systems different interpenetrating fluid components (eg. particles and entropy 
in heat conducting situations or a superfluid condensate and finite temperature excitations) can flow with distinct 
velocities. 

We have examined plane-wave propagation for a generic two-fluid system. By imposing the constraints of absolute 
stability (necessary in order for the components not to separate already in the absence of flow) and causality we 
have established limits on the equation of state (as represented by the master function A). In particular, we place 
constraints on the free sound speeds, the cross-constituent coupling and the entrainment. Some of the obtained 
results are more or less trivial extensions of the single fluid result, but others are unique to the two-fluid problem. 
The condition ()88|) on the entrainment is a new result, and as such serves as a new condition on, say, the kind of 
a ~ Lu model used by Comer and Joynt to model entrainment in the outer cores of neutron stars. 

We have demonstrated (for the first time) the existence of a relativistic two-stream instability. This is a generic 
phenomenon, that does not require particular fine-tuning to be triggered, nor is it limited to any specific physical 
system. The only requirement is that there is a relative (background) flow and some type of coupling between the 
fluids. While it is true that a single fluid can have an analogous instability, it is only active at an interface where 
there is shearing motion. Our analysis assumes that the two fluids are interpenetrating. 

In order to exhibit the generic nature of the two-stream instability, we have kept the analysis rather abstract. 
On the one hand, this means that it should be relatively straightforward to apply our results to particular physical 
systems. On the other hand, it means that we have not yet discussed the relevance of the instability for any particular 
system. A more detailed consideration of multi-fluid problems in relativity is required in order to establish whether 
this mechanism operates in nature. There are already exciting results that hint at this class of instabilities being 
associated with pulsar glitches [1, Q. In addition to exploring possible situations where these instabilities may 
operate, it would be very interesting to probe the nonlinear development of the unstable waves. So far, all studies 
have been at the linear perturbation level. The results establish the presence of the instability, but they do not 
shed any light on what happens once the unstable oscillation reaches nonlinear amplitudes. Detailed studies of this 
problem are essential if we are to understand the actual dynamical role of this instability. 
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